Two-Particle Bound States in the Continuum realized in the one-dimensional 

Hubbard Model with an Impurity 
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We report a bound state, which is embedded in the continuum spectrum, of the one-dimensional 
two-particle (Bose or Fermion) Hubbard model with an impurity potential. The state has the Bethe- 
ansatz form, although this model is nonintegrable. Moreover, for a wide region in parameter space, 
its energy is located in the continuum band. A remarkable advantage of this state with respect to 
similar states in other systems is the simple analytical form of the wave function and eigenvalue. 
This state can be tuned in and out of the continuum continuously. 
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A widespread preconception in quantum mechanics is 
that a normalizable (bound) state cannot be degener- 
ate in energy with a non-normalizable (extended) state. 
Generically, this is indeed the case: as argued by Mott 
in connection with the existence of sharp mobility edges 
[l[ , degeneracy between a localized and an extended state 
would be unstable against an infinitesimal perturbation 
converting the localized state into an extended one. How- 
ever, as pointed out by von Neumann and Wigner as early 
as in 1929 (2], the so-called bound state in the continuum 
(BIC), which is localized in space yet whose energy falls in 
the continuum band, does exist. Notably, von Neumann 
and Wigner took a constructive approach and found a 
simple, beautiful algorithm to construct such an exotic 
state together with the corresponding potential @, Q . 

It was almost forty years after von Neumann and 
Wigner's original work that the concept of BIC surfaced 
again. This time, Stillinger argued that BICs can exist 
in a two-electron atom Q , but not for physically relevant 
parameters. The first attempt to realize a BIC experi- 
mentally was taken by Capasso et al [5] , but their state 
was a defect mode in a minigap, not a true BIC. More re- 
cently, BICs do have been demonstrated experimentally 
in photonics 0, 0] using the analogy between optical sys- 
tems and quantum mechanics, and they are also theoret- 
ically predicted to exist in some other systems @- 13 1 . In 



spite of these progress, we note that generally BICs are 
fragile objects, and so far they have been found only in 
few systems. Some of them are simply protected b% 
metry [7] or rely on the separability of variables 
These BICs appear somewhat artificial. 

In this paper, we report the discovery of a BIC in 
a one-dimensional two-particle (Bose or Fermion) Hub- 
bard model with an impurity potential. As we shall see 
below, the model in question and the BIC found there 
have several desirable features in comparison with pre- 
vious models and the associated BICs. First, the model 
is simple yet nontrival, and most importantly, not ar- 
tificial, in contrast to many of the constructions above. 
Second, though the model is nonintegrable, half of the 



eigenstates can be shown to have the Bethe-ansatz form. 
They are distinguished from the diffractive states by a 
Z2-symmetry of the model: the Bethe (non-diffractive) 
states are odd under this symmetry, whereas the others 
are even. Some of the former are bound states, among 
them the BIC. This leads to explicit expressions of its 
wave function and eigenvalue, from which we see that 
the BIC can be tuned in and out of the continuum band 
by varying the model parameters. 

The model we investigate consists of two identical spin- 
less bosons (or two spin-^ fermions in the spin singlet 
space) in a one-dimensional infinite lattice. A defect is 
located at x = 0, leading to a local potential V. The 
two particles interact through a on-site interaction U. 
The wave function of the two particles is denoted as 
f(xi,X2), with — co < xi.2 < +oo being integers. The 
exchange symmetry requires f(xi,x 2 ) = f(x 2 ,x%) for 
bosons, which we shall consider in the following (the 
fermionic case is completely analogous). The Hamilto- 
nian is defined by its action on the wave function 

Hf(xi,x 2 ) = - ^ [f(xi + a, x 2 ) + f(xi,x 2 + a)} 

a=±l 

+ [V(S Xl , + S X2 , ) + U5 XUX2 ] f{xi,x 2 ).{l) 

Here V is the value of the impurity potential while U is 
the on-site interaction between the two particles, with the 
hopping strength set to unity. Here it should be stressed 
that in the absence of impurity potential or particle- 
particle interaction, the model is simply solvable. How- 
ever, in the presence of both impurity and interaction, 
it becomes nonintegrable even in the two-particle sector. 
By the unitary transform. f(x\, x 2 ) —> (— ) Xl+X2 f(xi, x 2 ), 
we can see that the eigenvalues and eigenstates of a 
Hamiltonian with parameters (V, U) are related to those 
of a Hamiltonian with parameters (— V, —U) in a simple 
way. Therefore, in the following we shall confine ourselves 
to the attractive case V < 0. 

We aim to solve the eigenvalues and eigenstates of this 
system, and especially the bound states, i.e., states in 
which both particles are localized in the vicinity of the 
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FIG. 1. The three continuum bands revealed by using the two 
quantities of D = [asi| + | £Pa | and Pj^t. The three horizontal 
lines are clearly visible in each panel. The eigenstates and 
eigenvalues E n are solved by exact diagonalization on a (2N + 
l)-site (N = 20) lattice and with open boundary condition. 
The values of the parameters are (V, U) = ( — 1.5, 5.5). 
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FIG. 2. The X1-X2 lattice decomposed into six regions by 
the :ri-axis, a?2-axis, and x\ = X2 line. Note that adjacent 
regions share the boundary between them. Especially, the 
origin belongs to all the six regions. 



impurity. To demonstrate that one of the bound states 
lies in the continuum, we have to obtain an overview 
of the extended states and identify the three continuum 
bands associated. This classification can be verified nu- 
merically using a finite size system (see Fig. [1}. The 
first band corresponds to two particles that are neither 
captured by the impurity nor bound together by the in- 
teraction between them. The impurity and interaction 
cause phase shifts but do not contribute to the energy 
of the state, which is just the sum of the kinetic energy 
of both particles, and thus this band covers the interval 
[—4, +4]. For the second band, one particle is captured 
by the impurity but the other is not. The energy of the 
first particle is —y/V 2 + 4 [l6[ while that of the other 
one lies between —2 and 2. This band covers therefore 

+ 2]. The third band corre- 
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sponds to a delocalized molecule state |17|. That is, the 
two particles are bound together by the interaction be- 
tween them and the composite moves as a whole on the 
lattice. The impurity causes phase shifts or local modifi- 
cations on the wave function but does not contribute to 
the energy. This band covers [— s/XJ 2 + 16, U] if U < 
or 



The presence of these three bands can be demonstrated 
numerically by using some quantities which reveal the 
distinct nature of the extended states. The first quan- 
tity is the sum of the distance of the two particles to 
the defect, D = \xi\ + \x2\- Suppose we choose a lat- 
tice of 2N + 1 sites with N sites on each side of the 
defect. For the first and third bands, since the par- 
ticles move through the lattice cither independently or 
bound together, (D) should be ~ N. For the second 
band, since always one particle is localized around the 
defect, (D) is expected to be ~ N/2. The other quantity 
is the probability of finding at least one particle out- 
side of a ball of radius R and centered at the defect. 

That is, P£a{f) = Ema X {|x 1 |,|x2|}>fll/( 2; 1^2)| 2 - Sup- 

pose R = N/2, it is easy to see that for the second and 
third bands, P^ t should be ~ 0.5, while for the first band 
it should be ~ 0.75. These predictions are readily veri- 
fied numerically. In Fig. [TJ we see how the three bands, 
although overlapping in energy, are separated by using D 
and P^ t . Moreover, their band edges coincide with the 
predicted values. The quantities D and P^ ut can also be 
employed to identify bound states (e.g., the BIC below): 
for bound states, the values of D and P^ ut should not 
increase with the size of the lattice, as long as N is large 
enough. 

The impurity potential and the interaction between 
the two particles are effective only on the three lines of 
xi = 0, x-2 = 0, and x\ — xi- Away from these lines, 
we have free particles. This observation suggests to try 
an ansatz for the eigenfunctions of Bethe-form which is 
characterized by just two parameters, k\ and &2- Specif- 
ically, in region Ii (0 < x\ < x%) in Fig. [2 the wave 
function is postulated to read 



f{xi,x 2 ) 



'2=0,1 

xexp[i(-) CT1 fcia;i +CTo + i(-) a2 k 2 X2-a } , (2) 



E E 

ITo=0,l tTi=0,l 0-2=0,1 



[U, y/U 2 + 16] if U > [17 1- 



where h(<To, o\, C2) = 4cro + 2o\ + 02 + 1. The wave 
function in regions I2 and I3 are defined similarly but 
with the A's replaced by B's and C's, respectively. The 
value of the wave function in the other three regions is 
determined by the condition f(x\,X2) — f(x2,Xi). In 
each region, we have eight different plane waves. The 
reason is that the interaction between the two particles 
can exchange their momenta, and the impurity potential 
can reverse the momentum of a particle. 

With the wave function in the form above, the eigen- 
value equation Hf = Ef, with E — — 2cosfci — 2cosfc2, 
is satisfied away from the three lines. Now we need 
to fulfill it also on the three lines. The linking con- 
ditions on the interfaces between adjacent regions give 
rise to a set of 24 homogeneous linear equations of the 
24 unknowns of A's, B's, and C's. The 24 x 24 coeffi- 
cient matrix, which depends on U, V, and fci,2, needs 
to be, and is indeed, singular to admit nontrivial solu- 
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tions. However, instead of dealing with the 24 x 24 ma- 
trix, we employ a simplification. Note that the system is 
reflection invariant with respect to the impurity. Defin- 
ing [Tf](xi, X2) — /(— xi, — X2), we have THT = H. 
Therefore, we can classify the eigenstates into even and 
odd ones with respect to the (parity) symmetry T. It 
turns out that the even-parity eigenstates (especially the 

I 



ground state) do not have the form of ©. However, the 
odd ones do. 

For the odd case, f(xi,x 2 ) = — /(— Xi, — X2), we need 
Bi = —Bg-i and At = — Cg-i. It is readily verified that 
the former implies the latter. In the end, the linking 
conditions reduce to a set of linear equations of Bi<i<±: 



(U + 2i(sin ki - sin k 2 )) \ Bi 
(U — 2i(sin k\ — sin k 2 )) I B4 



V 



(U + 2«(sinfci + sinfe 2 )) B 2 + 



(U - 2i(sin fci + sin fc 2 )) B 3 - 



2i sin ki 
V 

2i sin fei 



V 



2i sin fei 



V 



2i sinfci 



(B 1 +B 3 )J - (f/-2i(sinfci -sinfe 2 )) [B 4 + 
(B 2 + B 4 ^j - (U + 2i(sin k x -sin ha)) f-Bi - 
(B 2 + &t)J - (17 - 2 2 (sin h + sin k 2 )) (b 3 - 
(B1 + B3) J - (U + 2i(smh + sin k 2 j) [B 2 + 



V 



2i sin k 2 
V 

2i sin k 2 
V 

2i sin k 2 
V 

2i sin fee 



(B 3 + Bi) 
(B1+B2) 
(S 3 + B 4 ) 
■{B l +B 2 



= 0, (3a) 
= 0,(3b) 
= 0, (3c) 
= 0.(3d) 



The 4x4 coefficient matrix is always singular and has 
a two-dimensional nullspace. The reason for the value 
of two is the time-reversal symmetry of the Hamilto- 
nian. In terms of Eqs. ([2]) and (J3]), for fci.2 being 
real, if (Bi, B 2 , B 3 , B4) is a solution of then so is 
{B i ,B 3 ,B 2 ,B 1 ). 

Here some remarks are necessary. Suppose we allow 
both Bose and Fermi symmetry. Then Eq. (flj itself 
is invariant under both the exchange of x\ <-$■ x 2 and 
the reflection of xi t2 ~ > —xi t2 . The two symmetries di- 
vide the Hilbert space into four sectors. For the anti- 
symmetric (fermionic) sectors, the interaction is ineffec- 
tive and we have virtually free fermions in an impurity 
potential. The wave functions are in the Slater-form and 
hence also in the Bethe-form. Therefore, only in one of 
the four sectors, i.e. the symmetric (bosonic) sector with 
even parity, the wave functions have not the Bethe-form 
and are diffractive. The occurence of diffraction in a re- 
lated model defined on a continuous line was shown in 
the classic work by McGuire The interesting point 
here is that the symmetry T leads to a decomposition of 
the bosonic Hilbert space into two subspaces of which one 
shows no diffraction and can be therefore considered in- 
tegrable. This is confirmed by an analysis of the algebra 
of scattering matrices and the associated Yang-Baxter 
relations [I9j ]. 

We now proceed to study the odd-parity localized 
states. For localized states, it is easy to see that among 
the four coefficients -Bi<i<4, only one can be finite. The 
reason is the following: in region I 2 {x\ < < x 2 ), one 



of the two terms of e l 



and 



diverges to in- 



finity while the other converges to zero as \xj\ — > 00. 
Therefore, their prefactors cannot be finite simultane- 
ously Without loss of generality, suppose B\ = I while 
B 2 = B 3 = B4 = 0. The (unnormalized) wave function 



I 

in region I 2 then reduces to, 

f(x 1 ,x 2 ) = e lklXl+lk2X2 -e 



(4) 



Since in region I 2 , x\ goes to —00 while x 2 goes to +00, 
we need |e lfel | > I > |e ifc2 | for / to vanish at infinity. 
Note that we can also define — k 2 as k\, and —k\ as k 2 . 
To avoid counting a solution twice, we take the conven- 
tion that |e tfel+lfe2 | < 1. From Eq. (J3J), we get two inde- 
pendent equations for U, V, fci,2, 



(U + 2i(sin k\ - sin k 2 )) I 1 + 
(U + 2i(sin k x - sin fc 2 )) ( 1 - 



V 



2i sin fei 
V 



2i sin k 2 

From the linking conditions, we get furthermore 
V , V 



0, (5a) 
0. (5b) 



Ai = l 



A 5 = A 6 = 0, A 7 = 



2isinfci 2isinfci 
V 



2i sin k 2 



,A 8 = -1 



A 2 = A i = 0, (6a) 
V 



2i sin k 2 



(6b) 



In region Ix (0 < x% < X2), x\ and 2:21 = x 2 — x\ are 
independent variables taking values on Z + U0. We rewrite 
the wave function as 

f[X\ X 2 ) — j\ ie i ( kl+k2 ^ XlJrik2X21 -\- A 3e i ( k 2-k 1 )x 1 +ik 2 x 2 i 

Now for the wave function to be localized, we need 
|e lfel+lfc2 | < I and A 8 = 0. We thus get from j6b) 

V = 2ismk 2 = e lk2 -e~ ik2 . (8) 

Eq. ((Sb) is then satisfied. For (JS^,) to be satisfied, we 
need either V + 2i sin fej = or U + 2i (sin fei — sin k 2 ) = . 
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FIG. 3. In the region labeled with BIC, the state with energy 
Ebi [see is a bound state embedded in the [—4, +4] con- 
tinuum. In the region labeled with BOC, this state is a bound 
state below the [—4, +4] continuum band, and outside of all 
the three continuum bands. The boundary between the two 
regions is determined by the condition = —4. In the re- 
gion labeled with 2nd, the second odd-parity bound state with 
energy Et2 appears, which is outside the continuum bands. 
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FIG. 4. (Color online) Images of the squared wave func- 
tion of the bound state in the continuum, with (a) (V, U) = 
(-2,-0.5), E bl = -0.3284 and (b) (V, U) = (-2,-1.5), 
£?(,! = —0.7669. These states are obtained in a finite lat- 
tice by exact diagonalization, with 100 sites on each side of 
the impurity and open boundary condition. They agree with 
the analytic expressions in Eqs. (|4]), &), ((6b), and Q. 



But the former is not acceptable, otherwise the condition 
\e tkl | > 1 and | e 4,Cl +* fc2 1 < 1 can not be satisfied simulta- 
neously. Therefore, U + 2i(sinfei — smk 2 ) = 0, or 

V-U = 2isink 1 =e lkl -e~ ikl . (9) 

Defining Z\ = e lkl and z 2 — e lk2 , we have V = z 2 — z^ 1 , 

V — U = Z\ — z-j -1 , and E = —Z\ — z^ 1 — Z2 — z^ 1 , 
with 1 22 1 < 1 < \zi | < \z2 1 _1 . Instead of studying for 
what kind of (V, U) pairs, we have solutions of {z\,Z2) 
satisfying this condition, we do the inverse. Because V < 
by assumption, we have < z% < 1. Depending on the 
sign of Z\, we have two cases: 

(i) < z 2 < 1 < zi < z~ x . We get < V - U = 
Zi - zf 1 < z^ 1 - z 2 = -V. That is, 2V < U < 

V < 0. The energy of the wave function is — 
-\JV 2 + 4 - y/(V - U) 2 + 4. It is easy to prove that 
E b2 < {-4, - VV 2 + 4-2, - VU 2 + 16}. Consequently, 
this bound state is below all the three continuum bands 
and is thus not a BIC. 

(ii) -z^ 1 < z\ < -1 < < z 2 < 1. We get 

> V - U = zi - zf 1 > -z^ 1 + z 2 = V. (10) 

That is, V < U < 0. The eigenenergy of the wave func- 
tion is 

E bl = -y/V 2 + 4 + y/(V - U) 2 + 4. (11) 

Now, it is easy to show that > E bl > {U, -\JV 2 + 4 + 
2}, and thus E bl falls outside of the second and third 
band. But it can fall in the continuum band [—4, +4] to 
be an embedded eigenvalue. For example, for (V, U) — 
(-2,-0.5), E bl = -0.3284, which is inside the [-4, +4] 
continuum. The condition for this state to be a BIC is 
E b > —4. In Fig. [3j we have charted the regime where 
this condition is fulfilled. Note that this bound state 
exists whenever V < U < 0. However, only in a subset 



of this regime does it fall in the continuum. Its energy 
can be tuned continuously across the boundary —4. In 
Fig. [4] we have plotted its squared wave function for two 
sets of parameters. 

In conclusion, we have demonstrated the existence of 
a bound state in the continuum in a system of two inter- 
acting particles in an impurity potential. This discovery 
is actually a by-product of an investigation of the follow- 
ing well-motivated problem: for what value of U > will 
the (attractive) impurity potential no longer be able to 
bind the two bosons? It turns out that in this simple 
yet nontrivial problem a BIC appears, which is probably 
the simplest nontrival BIC as far as the wave function 
and eigenvalue are concerned. Moreover, it should be 
stressed that, unlike most BICs studied before which are 
generally one-body BICs, here we have a two-body BIC. 
It would be interesting to consider also three-particle or 
many-particle cases to see whether "partial integrabil- 
ity" persists and many-body BICs of simple form are 
possible. Finally, we would like to point out that in this 
system there exist other bound states (e.g. the ground 
state). However, they are all of even parity (diffractive) 
and located outside the continuum. 
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problem and stimulating discussions, M. Sekania and D. 
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ported in part by the Deutsche Forschungsgemeinschaft 
through TRR 80. 



[1] N. F. Mott, Adv. Phys. 16, 49 (1967). 
[2] J. von Neumann and E. Wigner, Phys. Z. 30, 465 (1929). 
[3] F. H. Stillinger, Phys. Rev. A 11, 446 (1975). 
[4] F. H. Stillinger, J. Chem. Phys. 45, 3623 (1966). 
[5] F. Capasso, C. Sirtori, J. Faist, D. L. Sivco, S-N. G. Chu, 
and A. Y. Cho, Nature (London) 358, 565 (1992). 



■5 



[6] D. C. Marinica, A. G. Borisov, and S. V. Shabanov, Phys. 

Rev. Lett. 100, 183902 (2008). 
[7] Y. Plotnik, O. Peleg, F. Dreisow, M. Heinrich, S. Nolte, 

A. Szameit, and M. Segev, Phys. Rev. Lett. 107, 183901 

(2011). 

[8] G. Ordonez, K. Na, and S. Kim, Phys. Rev. A 73, 022113 
(2006). 

[9] M. L. Ladron de Guevara and P. A. Orellana, Phys. Rev. 
B 73, 205303 (2006). 

[10] H. Nakamura, N. Hatano, S. Garmon, and T. Petrosky, 
Phys. Rev. Lett. 99, 210404 (2007). 

[11] N. Moiseyev, Phys. Rev. Lett. 102, 167404 (2009). 

[12] J. W. Gonzalez, M. Pacheco, L. Rosales, and P. A. Orel- 
lana, Europhys. Lett. 91, 66001 (2010). 



[13] M. I. Molinal, A. E. Miroshnichenko, and Y. S. Kivshar, 

Phys. Rev. Lett. 108, 070401 (2012). 
[14] A. K. Jain and C. S. Shastry, Phys. Rev. A 12, 2237 

(1975). 

[15] M. Robnik, J. Phys. A 19, 3845 (1986). 

[16] R. P. Feynman, R. B. Leighton, and M. Sands, The 

Feynman Lectures on Physics Volume 3, Addison Wesley 

Longman (1970). 
[17] K. Winkler, G. Thalhammer, F. Lang, R. Grimm, J. 

Hecker Denschlag, A. J. Daley, A. Kantian, H. P. Biichler, 

and P. Zoller, Nature (London) 441, 853 (2006). 
[18] J. B. McGuire, J. Math. Phys. 5, 622 (1964). 
[19] J. M. Zhang, D. Braak, and M. Kollar (in preparation). 



